First glances at the transversity parton distribution through 
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We present first observations of the transversity parton distribution based on an analysis of 
pion-pair production in deep inelastic scattering off transversely polarized targets. The extraction 
of transversity relies on the knowledge of dihadron fragmentation functions, which we take from 
electron-positron annihilation measurements. This is the first attempt to determine the transversity 
distribution in the framework of coUinear factorization. 

PACS numbers: 13.88. +e, 13.87.Fh, 13.60.Hb 



The distribution of quarks and gluons inside hadrons 
can be described by means of parton distribution func- 
tions (PDFs). In a parton-model picture, PDFs describe 
combinations of number densities of quarks and gluons 
in a fast-moving hadron. The knowledge of PDFs is cru- 
cial for our understanding of Quantum Chromodynamics 
(QCD) and for the interpretation of high-energy experi- 
ments involving hadrons. 

If parton transverse momentum is integrated over, in 
the Bjorken Umit the partonic structure of the nucleon is 
described in terms of only three PDFs: the well known 
unpolarized, /f (x), and heUcity, gfix), distribution func- 
tions, and the transversity distribution function 
which measures the transverse polarization of quarks 
with flavor q and fractional momentum a; in a trans- 
versely polarized nucleon Intuitively, helicity and 
transversity give two orthogonal pictures of the partonic 
structure of polarized nucleons. They have very different 
properties, and transversity is much less known. In this 
work, we present an extraction of the hf. 

Transversity is related to the interference of amplitudes 
with different helicities of partons and of the parent nu- 
cleon. In jargon, it is called a chiral-odd function. There 
is no transversity for gluons in a nucleon, and has a 
pure non-singlet scale evolution Si]. From transversity 
one can build the nucleon tensor charge, which is odd 
under charge conjugation and can be computed in lattice 
QCD [3] (for a review on transversity, see Ref. [A] and 
references therein). 

Transversity is particularly difficult to measure be- 
cause it must appear in cross sections combined with 
another chiral-odd function. An example is the cross 
section for single-particle inclusive Deep Inelastic Scat- 
tering (DIS), where hi appears in a convolution with 
the chiral-odd Collins fragmentation function iJj^' [5], 
which describes the correlation between the transverse 
polarization of a fragmenting quark with flavor q and the 
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transverse momentum distribution of the detected un- 
polarized hadron. The convolution hi (g) gives rise 
to a specific azimuthal modulation of the cross section. 
The amplitude of the modulation has been measured by 
the HERMES and COMPASS collaborations [5]. In or- 
der to extract the transversity distribution from this sig- 
nal, the Collins function should be determined through 
the measurement of azimuthal asymmetries in the dis- 
tribution of two almost back-to-back hadrons in e^e~ 
annihilation [7|. The Belle collaboration has measured 
this asymmetry making the first-ever extraction of 
h\ possible from a simultaneous analysis of ep^ e'wX 
and e+e^ ttttA' data [H]. 

In spite of this achievement, some questions still hin- 
der the extraction of transversity from single-particle- 
inclusive measurements. The most crucial issue is the 
treatment of evolution effects, since the measurements 
were performed at very different energies. The convo- 
lution hi (E) involves the transverse momentum of 
quarks. Hence, its evolution should be described in the 
framework of the transverse-momentum-dependent fac- 
torization |10| . Quantitative explorations in this direc- 
tion suggest that neglecting evolution effects could lead 
to overestimating transversity [llj . 

In this context, it is of paramount importance to ex- 
tract transversity in an independent way, requiring only 
standard coUinear factorization where the above com- 
plications are absent (see, e.g, Refs. [T2] and references 
therein). Here, we come for the first time to this result by 
considering the semi-inclusive deep- inelastic production 
of two hadrons with small invariant mass. 

In this case, the transversity distribution function 
is combined with a chiral-odd Dihadron Fragmentation 
Function (DiFF), denoted as H^'^ [13], which describes 
the correlation between the transverse polarization of the 
fragmenting quark with fiavor q and the azimuthal ori- 
entation of the plane containing the momenta of the de- 
tected hadron pair. Contrary to the Collins mechanism, 
this effect survives after integration over quark trans- 
verse momenta and can be analyzed in the framework 
of coUinear factorization. This process has been studied 
from different perspectives in a number of papers p!3Hl6] . 
The only published measurement of the relevant asymme- 
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TABLE I: Semi-inclusive DIS data of the asymmetry ylois 
from HERMES 17^. The errors are mainly statistical (we 
added the systematic errors in quadrature). The average val- 
ues of the variables are taken from Tab. 5.1 of Ref. [23]. The 
other variables have been integrated in the range 0.5 < Mh < 
1 GeV and 0.2 < ^ < 1. 



try has been presented by the HERMES collaboration for 
the production of tt+tt" pairs on transversely polarized 
protons |17| . Preliminary measurements have been pre- 
sented by the COMPASS collaboration [H]. Related pre- 
liminary measurements in proton-proton scattering have 
been presented by the PHENIX collaboration P^ . 

Similarly to the single-hadron case, we need to in- 
dependently determine H"^ ^ by looking at correlations 
between the azimuthal orientations of two pion pairs in 
back-to-back jets in e+e" annihilation The mea- 

surement of this so-called Artru-Collins azimuthal asym- 
metry has recently become possible thanks to the Belle 
collaboration [22* . 

In the HERMES publication flT], the asymmetry was 
denoted as sine. brevity and without am- 

biguity, here we will use the notation v4dis- The data 
set was collected in bins of the variables x (the momen- 
tum fraction of the initial quark), z (the fractional energy 
carried by the tt+tt" pair), and Mh (the invariant mass 
of the pair). Since our interest here lies mainly on the 
transversity distribution, and to avoid problems when 
dealing with three different projections of the same data 
set, we consider only the x binning. In Tab. |Ij we repro- 
duce the data for convenience indicating for each bin also 
the average hard scale (Q^) and fractional beam energy 
loss {y). 

The ^Idis measured by HERMES [T7] can be inter- 
preted as P3] 



Amsix^Q^) = ~Cy 



E,e2/i?(x,Q2)„T(Q2) 



E,e2A«(x,Q2)„^(Q2) 

where (neglecting target-mass corrections) 

{1-y) ^ l-{y) 
^ (l-y + yV2) ~ l-(y> + (2/>V2 ' 

In Eq. ([I]), we also introduced the following quantities 
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with \R\/Mh = ^1/4- ml/M^ . Dl'^'' is the unpo- 
larized DiFF describing the hadronization of a quark q 
into a 'k^t:~ pair plus any number of undetected hadrons, 
averaged over quark polarization and pair orientation. 

Finally, H^l^^ ^ is a chiral-odd DiFF, and denotes the 

component of H^'^^^ ^ that is sensitive to the inter- 
ference between the fragmentation amplitudes into pion 
pairs in relative s wave and in relative p wave, from which 
comes the common name of Interference Fragmentation 
Functions [15 . Intuitively, if the fragmenting quark is 
moving along the z direction and is polarized along y, a 

positive H^l^^ means that tt"*" is preferentially emit- 
ted along ~x and tt^ along x. Since in this case no 
ambiguities arise, in the following we shall conveniently 
simplify the notation by using Df and H^'^ to denote the 
relevant DiFFs. 

In our analysis, we make the following assumptions 
(valid only for tt+tt" pairs) based on isospin symmetry 
and charge conjugation |24) : 

D'l = Df = Dl = Df , (4) 

Dl^Dl, Dl^Dl, (5) 

Hf"" = -Hf^ = -Hf" = Hf^ , (6) 

= -Hf ^ Hf' = -Hf"" = . (7) 

We also assume Dl = Ns -D" and we consider the two 
scenarios A^s = 1 and A^s — 1/2. The second choice is 
suggested by the output of the PYTHIA event genera- 
tor [35]. Our final results will not depend strongly on 
this choice. 

The above assumptions allow us to turn Eq. ([l]) into 
the following simple relation (neglecting charm quarks) 



xhl'^{x,Q^ 



\xh'i-{x,Q^) = - 

nl(g2 



^Dis(a:,Q^) 

Cy 

e^N (8) 



sp 



{z,MtQ'), (3) 



where A^„ = A^^ = 1 and f^+" = ff + /f , ^h\- h\. 
Our goal is to derive from data the difference between 
the valence up and down transversity distributions by 
computing the r.h.s. of the above relation. 

The PDFs in Eq. ^ can be estimated using any 
parametrization of the unpolarized distributions. We 
chose to employ the MSTW08LO PDF set [26 . We 
checked that using different sets makes no significant 
change. We also checked that the charm contribution 
is irrelevant. 

The only other unknown term on the r.h.s. of Eq. ([S]) 
is the ratio n„/nj^. We extract this information from 
the recent measurement by the Belle collaboration (22] of 
the Artru-Collins azimuthal asymmetry ^cosC^a-k^b) ^ 
1211 [23 (denoted as ai2R. in the experimental paper). As 
in the previous case, without ambiguity we simplify the 
notation and refer to this asymmetry as A^+e-- Using 
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the assumptions Q-Q, the asymmetry can be written 



A,+,^{z,Ml,z,Ml,Q^) 



(sm( 



(sin( 



M 



Ah 



(1 + cos2 02) (5 + 7V|) DUz, MlQ^) D^iz, Ml Q^) +4D-{z, M^Q^) D1{z, Ml Q^) 



r 



where 6*2 is the angle between the direction of the lepton 
annihilation and the thrust axis, and 9 is the angle be- 
tween the momentum of one hadron in the cm. of the 
hadron pair and the total momentum of the pair in the 
laboratory JTB] . 

In Eq. (l8|) we need n^/uu at the experimental val- 
ues of (Q^) of Tab. |l] and integrated over the HERMES 
invariant-mass range 0.5 < M^ < 1 GeV. We will get 
to this number in two steps: first, we estimate the ratio 
n\/nu integrated over 0.5 < Mh < 1 GeV and at the 
Belle scale (100 GeV^), then we address the problem of 
changing Q^. 

We consider the Belle asymmetry integrated over (z, z) 
and binned in {Mh,Mh)- We restrict our attention only 
on the bins between 0.5 < {Mh,Mh) < 1-1 GeV (neglect- 
ing the small difference with the HERMES upper limit). 
We weight the contribution of each bin by the inverse of 
the statistical error squared, which should be to a good 
approximation proportional to the denominator of the 
asymmetry in each bin. By summing over all bins in the 
considered range, we get the total asymmetry 



A, 



(sin^ 02) 



(sing) (sing) 5(0^ 

(l + cos2 6*2) (5 + 7V|)n2 +4n2 
-0.0307 ±0.0011 . 



(10) 



From the Belle analysis, we know that 
(sin2 02) 



(1 + COs2 



0.753 , (sing) (sing) = 0.871 



4n2 



(5 + iV2)n2 
Therefore we obtain 



= 0.415 ±0.047 



4/n„(100 GeV^) ^ -0.273 ± 0.007cx ± 0.009th , (H) 

where the second error comes from using the two different 
values of the s— quark normalization N^. We assumed 
the sign of the ratio to be negative in order to obtain a 
positive M— quark transversity distribution. To verify the 
reliability of this procedure, we repeated the calculation 
estimating the denominator of the asymmetry using the 
PYTHIA event generator |25j without acceptance cuts 
(courtesy of the Belle collaboration). The result falls 
within the errors quoted in Eq. ( 11 ). 

The last step in the procedure is to address the 
(52— evolution of nj^/n^, since the Belle scale is very dif- 
ferent from the HERMES one (see Tab. |l|. DiFFs must 



be connected from one scale to the other via their QCD 
evolution equations [55]. In order to do this, we need 
to know the z dependence of H^'^ and for each Mh 
value. For Hf^, we fit the Belle data for A^+e- binned 
in (z, Mh) and integrated over (z, M^), multiplied by the 
inverse of the statistical error squared. 

The D\ should be obtained from global fits of unpolar- 
ized cross sections, similarly to what is done for single- 
hadron fragmentation functions [29]. In the absence of 
published data, we extract D\ by fitting the unpolarizcd 
cross section as produced by the PYTHIA event gener- 
ator [3S], which is known to give a good description of 
the total cross section. Following the assumptions intro- 
duced in Eqs. Q and ([s]), we describe the unpolarized 
cross-section for the production of a hadron pair with [21] 



da 



47ra2 



dzdMl Q2 



9^9^ 



9 ^ 



(12) 



where a is the fine structure constant. We assume the 
integration over cos g2 to be complete in the Monte Carlo 
sample. 

We start from a parametrization of D\ and Hf^ at 
Qq = 1 GeV2. Taking inspiration from the model analy- 
sis of Ref. fM' , for both DiFFs we consider two channels 
which are effective in the considered range 0.5 < Mh < 
1.1 GeV: the fragmentation into the p resonance decaying 
into TT+TT", and the continuum arising from the fragmen- 
tation into an incoherent tt+tt" pair. Then, we evolve 
the DiFFs at LO using the HOPPET code [SD], which 
we suitably extended to include also chiral-odd splitting 
functions. Finally, we fit the cross section (12 1 and the 
numerator of the asymmetry ([9]) in the bins of interest. 
We checked that the final results are affected in a neg- 
ligible way by the gluonic component Df{z,Mh;Ql). A 
thorough analysis will be presented in a future publica- 
tion 

By integrating the extracted DiFFs in the HERMES 
range 0.5 < Mh < 1 GeV and 0.2 < z < 1, we can 
calculate the evolution effects on n\/nu at each (Q^) in- 
dicated in Tab. HJ It turns out that the ratio is decreased 
by a factor 0.92±0.08, where the error takes into account 
the difference of in the HERMES experimental bins 
as well as the uncertainty related to different starting 
parametrizations at Qg = 1 GeV^. In conclusion, for the 
extraction of transversity in Eq. ([S]) we use the number 
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FIG. 1: The ar/i" " - a;/if" /4 of Eq. Q as a function of x. The 
error bars are obtained by propagating the statistical errors of 
each term in the equation. The uncertainty band represents 
the same observable as deduced from the parametrization of 
Ref. El. 



In Fig. |T] the datapoints denote the combination 
xft."^ - of Eq. (jS), plotted for each {x) and (Q^) 

listed in Tab. lU We studied the influence of the errors of 
each element in the r.h.s. of Eq. ([s]). The only relevant 
contributions come from the experimental errors in the 
measurement of ^dis, as reported in Tab. |lj and from 
the 9% theoretical uncertainty on nu/n\. 

In Fig. [T] the central line represents the best fit for the 
combination x/i"" ~ x/i^"/4, as deduced from the most 



recent parametrization of /i"" and /ij^" extracted from the 
Collins effect [32]. The uncertainty band is obtained by 
considering the errors on the parametrization and tak- 
ing the upper and lower limits for the combination of 
interest. Our data points seem not in disagreement with 
the extraction. However, a word of caution is needed 
here: while the error bars of our data points correspond 
to Ifj deviation from the central value, the uncertainty 
on the parametrization [32j corresponds to a deviation 
Ax^ ~ 17 from the best fit (see Ref. [33] for more de- 
tails) . In any case, to draw clearer conclusions more data 
are needed (e.g., from the COMPASS collaboration 18!). 

In summary, we have presented for the first time a 
determination of the transversity parton distribution in 
the framework of collinear factorization by using data for 
pion-pair production in deep inelastic scattering off trans- 
versely polarized targets, combined with data of e^e~ 
annihilations into pion pairs. The final trend of the ex- 
tracted transversity seems not to be in disagreement with 
the transversity extracted from the Collins effect |32j . 
More data are needed to clarify the issue. 
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